Strong Effects of Cation Vacancies on the Electronic and Dynamical Properties of FeO 
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We report pronounced modifications of electronic and vibrational properties induced in FeO by 
cation vacancies, obtained within density functional theory incorporating strong local Coulomb in- 
teractions at Fe atoms. The insulating gap of FeO is reduced by about 50% due to unoccupied 
electronic bands introduced by trivalent Fe ions stabilized by cation vacancies. The changes in the 
electronic structure along with atomic displacements induced by cation vacancies affect strongly 
phonon dispersions via modified force constants, including those at atoms beyond nearest neighbors 
of defects. We demonstrate that theoretical phonon dispersions and their densities of states repro- 
duce the results of inelastic neutron and nuclear resonant x-ray scattering experiments only when 
Fe vacancies and Coulomb interaction U are both included explicitly in ab initio simulations, which 
also suggests that the electron-phonon coupling in FeO is strong. 

PACS numbers: 71.30.-(-h, 63.20.dk, 63.20.kp, 71.23.An 



As explained by the theory of Motti and Hubbard;^ the 
insulating state in transition metal oxides results from 
the electron localization in the narrow 3d band due to the 
local Coulomb interactions Iron oxides are in the focus 
of present-day research covering a broad range of complex 
phenomena such as the Verwey transition in magnetite^ 
or structural instabilities in polar nanosystemsi^ High- 
resolution neutron and x-ray scattering experiments on 
magnetite^ suggest a considerably more complex charge 
order than originally invoked by Verweyi^ It has been 
shown that the microscopic understanding of this phe- 
nomenon requires full information about both the elec- 
tronic structure and the lattice dynamics in the presence 
of strong electron interactions which significantly amplify 
the electron-phonon coupling, as predicted in the theory^ 
and confirmed recently by inelastic x-ray scatteringi^ 
This could be a common feature of iron oxides — in- 
deed, the exchange-induced phonon splitting observed in 
wiistitei^ suggests a strong coupling between electronic 
and lattice degrees of freedom. 

Wiistite, Fei_a;0, is another strongly correlated insu- 
lator within iron oxides, with the insulating energy gap of 
~ 2.4 eV arising mainly from the 0(2p) — > Fe(3(i) charge- 
transfer processiii Below the Neel temperature T/v = 198 
K antiferromagnetic (AF) order sets in, which suggests 
that also here strong on-site Coulomb interaction U plays 
a role. A stable crystal of Fei-^O shows large Fc- 
deficiency on the level 0.05 < a; < 0.15^^ which substan- 
tially influences its structural properties] ^^d^ It may also 
affect the geophysical processes in the Earth's lower man- 
tle where wiistite, existing in the (Mg,Fe)0 solid solution, 
is one of the most important constituents So far, first- 
principles investigations on a vacancy-defected FeCU^ 



and other transition metal oxidesi^i^ are scarce. Previ- 
ous studies were devoted to the stoichiometric FeO )^^i^° 
leaving the effect of cation vacancies on the electronic 
and lattice properties of wiistite largely unexplored. 

At room temperature, Fei_2;0 crystallizes in a rock- 
salt structure and exhibits the rhombohedral distortion 
with a weak elongation along the [111] direction in the 
AF phase below T^. The magnetic moments on the Fe 
cations are aligned perpendicular to the ferromagncti- 
cally ordered (111) planes reversing their orientations 
from one to another planei^i Due to the cation defi- 
ciency and crystal-field effects, the average magnetic mo- 
ment m = measured at 4.2 K is significantly 
reduced from the ionic Fe^+ value m = AjiB (for spin 
5 = 2)<2^ Large splitting of phonon modes below Tm 
indicates strong spin-phonon coupling in wiistite 

In this Rapid Communication we demonstrate that 
Fe vacancies (Vpe) strongly modify the structural, elec- 
tronic, and vibrational properties of wiistite. Our the- 
oretical investigations uncover qualitative difference be- 
tween stoichiometric FeO and defected Fei-^^O contain- 
ing cither 3% or 6% of cation vacancies. The present 
ah initio results for Fei_a;0 are in good agreement with 
the available experimental data from inelastic neutron 
scattering (INS)^ and nuclear resonant inelastic x-ray 
scattering (NRIXS)^^ which validates our findings. 

Notably, a realistic description of the electronic struc- 
ture of Mott insulating states in transition metal monox- 
ides requires methods which implement explicitly lo- 
cal Coulomb interactions Previous theoretical stud- 
ies within the local density approximation (LDA) using 
the Coulomb interaction U (LDA-|-C/)^i2^ confirmed the 
prominent role of local electron interactions on the band 



2 



structure of defect-free FeO and significant contribution 
of the 0{2p) states at energies just below the Fermi en- 
ergy. Cation vacancies are responsible for charge redis- 
tribution and appearance of trivalent Fe ions in the oc- 
tahedral and tetrahedral (interstitial) positions J^i^^ The 
highly correlated nature of 3d electronic states has also a 
strong impact on the lattice dynamics of transition metal 
oxides. The experimental data have been well reproduced 
within the LBA+U approach^^ applied to NiO^ CoO^ 
and MnOj^ further improvements were obtained in the 
Green's function approachi^ Until now, lattice dynamics 
of stoichiometric FeO has only been studied in a semiem- 
pirical model)^ 

The present calculations were performed within the 
generalized gradient approximation (GGA) and the pro- 
jector augmented-wave method^l implemented in the 
VASP code<^ The on-site interactions at Fe ions are: 
the Coulomb element U = 6 eV and Hund's exchange 
J ~ 1 cVj^i21 The AF superccU composed of 64 atoms 
was used to calculate stoichiometric FeO. The Fei_a;0 
crystals with x = 3% and x = 6% were modeled by re- 
moving either one or two Fe atoms from the optimized 
AF supercell. Various Vpc-Vpe distances within the sim- 
ulated supercell with x = 6% were considered. Finally, 
a configuration with the lowest energy, corresponding to 
the distance between two vacancies of 1.22a, with a be- 
ing the lattice constant of the perfect FeO, was chosen 
for further band structure and phonon calculations. The 
energy difference between the lowest and the next stable 
configuration amounts to 7 meV/atom indicating rather 
high mobility (diffusion) of vacancies in wiistite. The vol- 
umes of the supercells representing Fei_a;0 systems were 
fixed at the optimized value of the vacancy-free lattice, 
the symmetry constraints were removed, and all atoms 
were allowed to relax. 

We identified two distinct valence states of Fe cations 
for each Fei_a.O composition. Majority of cations are 
Fe2+ ions with moments m = 3.72/xb, but there are two 
iron ions per each Vpe that have enhanced charges and 
moments (m = 4.22/1^). For single- vacancy system, they 
are located at distances a and 1.73a from Vpe- The Bader 
charge analysis^^ indicates an increase in their valence 
charges by -|-0.31e over the valence of Fe^+ ions. We 
assign the above features to the Fe'^"'" ions. 

The orbital projected electron densities of states 
(DOS) for FeO and the two Fci-^-O systems considered 
here arc dominated by the 0(2p) and Fc{3d) states; the 
latter showing the crystal-field splitting between the t2g 
and Cg symmetry. As in the previous studieSf^^ the elec- 
tronic states close to the energy gap are occupied by the 
minority spin t2g electrons, see Fig. [TJ The empty Fe(3(i) 
states are split and shifted to higher energies due to the 
Hubbard interaction U, while the 0{2p) states are almost 
fully occupied and hybridize with the Fe(3(i) states in a 
broad energy range below the Fermi energy. 

We observe significant modifications of the electronic 
structure for the system with incorporated cation vacan- 
cies. First of all, the projected electronic DOS becomes 
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FIG. 1. (Color online) Total orbital projected electronic 
DOS obtained for: (a) the stoichiometric FeO, and defected 
Fei_a;0, with (b) x = 3%, and (c) x = 6%. The positive (neg- 
ative) DOS represents the spin-up (spin-down) states in each 
case; the gaps between the valence and conduction band are 
indicated by dashed lines. Parameters: 17 = 6 eV, J = 1 eV. 



asymmetric in the distribution of the spin-up and spin- 
down components of the 3d states after removing one Fe 
cation from the spin-down sublattice. Also, unlike in the 
stoichiometric FeO, see Fig.[T](a), the oxygen bands of the 
wiistite crystals with different concentration of Fe vacan- 
cies remain spin-polarized. Each Vpe induces changes in 
the charge distribution of its immediate neighborhood, 
pushing the 0(2p) states closer to the Fermi energy. As 
the most interesting feature we identify additional elec- 
tronic bands that originate from the 3d states of Fe'^+ 
ions. In Fei_a;0 with x = 3%, the occupied 3d states 
of Fe'^+ ions form a narrow band at about —7.6 eV, see 
Fig. [IKb). It lies below the main manifold of the Fe^"*" 
states and consists of all five Fe(3d) states. The unoccu- 
pied Fe'^"'"(3d) states appear as the lowest energy empty 
states and reduce the insulating energy gap of wiistite 
from ~ 2.2 eV to ^ 1.8 eVM- Spin polarization of the 
Fe'^^ states depends on the Fe (spin-up or spin-down) 
sublattice into which the Vpo is introduced. The elec- 
tronic bands due to the trivalent Fe ions in the Fei_a;0 
system with x = 6% exist in both spin channels as the 
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FIG. 2. (Color online) Phonon dispersion relations for sto- 
ichiometric FeO (solid lines) and the relative intensities of 
the phonon modes in Fei-^O with x — 3%, for U = 6 eV, 
J = 1 eV. The experimental data (solid squares) are adopted 
from Ref.— . The highest intensity mode is taken as a ref- 
erence (100%); relative intensities are arranged into the fol- 
lowing ranges indicated by the intensity of gray scale (color): 
5-20% ( , 20-50% , 50-80% (), 80-100% (). High-symmetry 
points are (in units of 2-K/a): F — (0,0,0), X = (^,^,0), 
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cation vacancies occupy both spin-sublattices. A larger 
concentration of vacancies x =^ 6% broadens electronic 
bands and leads to further decrease of the insulating gap 
to the value of 1.3 eV— . This value is quite close to the 
observed optical gap of 1.0 eV.— The electronic states 
induced by disorder are expected to influence transport 
properties of FeO<^ 

Phonons in FeO and Fei_2;0 were determined by the 
direct method,— vifhere the Hellmann-Feynman forces 
were obtained by displacing all symmetry non-equivalent 
atoms from their equilibrium positions with the ampli- 
tude of 0.03 A and applying both positive and nega- 
tive displacements to minimize systematic error <^ The 
transverse optic (TO) phonon modes followed directly 
from the diagonalization of the dynamical matrix)^ while 
the longitudinal optic (LO) modes were obtained by in- 
troducing into the dynamical matrix the non-analytical 
term;^ the latter dependent on the Born effective charge 
tensor Z* and the electronic part of dielectric constant 
Eoo- Both Z* and eoc were determined within the lin- 
ear response method^ implemented in the VASP code. 



Our calculations resulted in 



= 2.2 and 



Eoo = 5.2624^ Ctxj and \Z*\ averaged over particular 
sublattices remain almost insensitive to the incorporated 
cation vacancies?^ whereas the Fe"^"*" ions acquire the ef- 
fective dynamical charge of -1-1.3 over that for Fc2+ ions. 

Imposed cubic symmetry constraints result in 6 
phonon branches. The transverse acoustic (TA) and TO 
phonon modes remain doubly degenerate along the F- 
X and T-L directions. Despite an apparent correlation 
between our theoretical curves and the INS data^ we 
notice some discrepancies for the TO branch close to the 
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FIG. 3. (Color online) Partial phonon DOS obtained for stoi- 
chiometric FeO (solid lines) and Fei_a;0 with x = 6% (shaded 
area) for sublattice of: (a) oxygen and (b) iron ions. Experi- 
mental points in (b) are taken from the NRIXS experiment.— 
Dotted lines denote the Gaussian convolutions of our theo- 
retical results with the experimental resolution of 2.4 meV. 
Parameters: U — 6 eV, J = 1 eV. 



L point and the longitudinal acoustic (LA) branch near 
the X point, see Fig. [5] An evident departure of the cal- 
culated LO branch from the experimental points along 
the T — X direction is mainly due to the interpolation 
function used here to calculate the LO-TO splitting42. 

Some of the above discrepancies can be explained by 
considering phonons in the FeO crystals defected by 
cation vacancies. However, we have to remark that an 
elaboration of the resulting phonon dispersion curves in 
the system with point defects requires a special treatment 
which involves application of the phonon form factor4^ 
We follow the same methodology as used before for 
the strongly correlated CoO with various point defects 
(cation vacancies and Fe impurities) The resulting in- 
tensities of the phonon modes in Fei_a;0 with x = 3% 
are overlaid onto the phonon dispersions for the stoichio- 
metric FeO in Fig. [21 Discontinuities in some phonon 
branches come from the 5% cut-off applied to the relative 
intensities of the phonon modes in Fei_a;0. Although the 
calculated phonon spectrum is dominated by the modes 
of low relative intensities, the INS data concentrate close 
to the modes with intensities exceeding 30%. Moreover, 
the calculated phonons at the Brillouin zone boundaries, 
e.g. at the L and X points, match better the experi- 
mental data. Finally, despite many similarities between 
the phonon spectra of FeO with 3% and 6% vacancies, 
the system with doubled vacancy concentration exhibits 
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FIG. 4. (Color online) Diagonal elements $ of the on-site 
force-constants matrix for distance d/a from the vacancy in 
the Fei_a;0 system with x = 3%. Vertical dashed lines denote 
ideal positions of the Fe and O atoms, whereas the horizontal 
solid and dashed lines stand for the respective force constants 
in the stoichiometric FeO. Parameters: (7 = 6 eV, J — 1 eV. 



significantly more phonon modes carrying low intensities. 

Essential differences in the dynamics of the stoichio- 
metric and defected lattices are revealed by the phonon 
DOS, see Fig. [31 Here, our theoretical results of par- 
tial Fe DOS obtained for stoichiometric FeO and Fei_2;0 
with a; = 6% are compared to the experimental data of 
NRIXS^ performed on the wiistite sample with the non- 
stoichiometry of 5.3%, see Fig.[3Jb). The iron and oxygen 
contributions span different energy ranges which over- 
lap weakly due to large difference between their masses. 
Both convoluted and bare theoretical phonon spectra of 
stoichiometric FeO show quite narrow and well-resolved 
peaks, opposite to the broad phonon maxima in the 
Fei_a:0 spectra. These latter spectra agree remarkably 
well with the experimental results A smearing of the 
two-peak structure, inevitably connected with the pres- 
ence of Fe vacancies, becomes even more enhanced for 
the 0-sublattice, indicating its greater sensitivity to the 
changes in the local crystal environment. 

To clarify the origin of changes in atomic vibrations, we 
have analyzed the force constant matrix for the wiistite 
crystal with x = 3%, see Fig. 01 The on-site force con- 
stants <I>a along the main crystallographic axes a = x^y, z 
increase for the vacancy neighboring oxygen atoms over 



the respective ones in the idealized vacancy- free FeO. 
They are displaced outward the vacancy by 0.18 A. On 
the other hand, the O atoms occupying the subsequent 
shells are attracted toward the nearest Fe atoms and their 
force constants remain equally distributed around the 
value corresponding to the perfect lattice. Similar but 
weaker effects on force constants have also been encoun- 
tered in the cation-deficient CoOJ^ 

Relaxations of the atomic positions in the Fe-sublattice 
are less pronounced, except the second nearest neighbors 
to a vacancy which shift closer to it4i The force constants 
at the iron sites show considerably broader distribution 
range than those at oxygens sites. Moreover, the force 
constants at Fc'^"'' ions are significantly enhanced from 
their values at Fe^+ ions, with the largest values of 
at the distance a from Vpo- This clearly demonstrates a 
strong impact of the wiistite electronic structure on its 
interatomic forces and accounts for the observed broad- 
ening of peaks in the phonon DOS. It may also explain 
large spectral linewidths of the dielectric loss function 
and strong electron-phonon coupling observed recently 
for FeO with the nonstoichiomctry of 7-8% Such pro- 
nounced changes in the phonon DOS, in particular the 
excess of spectral intensity at lowest energies, infiuence 
also thermodynamic properties of FeO. For instance, in 
case of one cation vacancy the calculations predict the 
increase of lattice heat capacity at T = 50 K by about 
10% over the stoichiometric FeO, and enhanced thermal 
atomic displacements in the neighborhood of a vacancy. 

Summarizing, we have reported dramatic changes of 
the electronic structure and vibrational spectra of wiistite 
caused by cation vacancies. Significant charge redistribu- 
tion on Fe ions induces the empty electronic states above 
the valence band and strongly reduces the insulating 
gap. The changes in the electronic structure and atomic 
displacements generated by defects substantially modify 
the force constants not only in the immediate neighbor- 
hood of vacancies but also at more distant atoms. It 
results in broadening of peaks in the phonon spectra, in 
close agreement with available experimental data. Pre- 
dicted changes in the dynamics of the oxygen sublattice 
in Fei_2;0 provide an experimental challenge that could 
be resolved, for example, by incoherent INS experiments. 
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